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We investigate the structure of halos in the sDGP (self-accelerating branch of the Dvali- 
Gavadadzc-Porrati brancworld gravity) model and the galilcon modified gravity model on 
the basis of the static and spherically symmetric solutions of the collisionlcss Boltzmann 
equation, which reduce to the singular isothermal sphere model and the King model in 
the limit of Newtonian gravity. The common feature of these halos is that the density of 
a halo in the outer region is larger (smaller) in the sDGP (galileon) model, respectively, 
in comparison with Newtonian gravity. This comes from the suppression (enhancement) 
of the effective gravity at large distance in the sDGP (galileon) model, respectively. 
However, the difference between these modified gravity models and Newtonian gravity 
only appears outside the halo due to the Vainshtein mechanism, which makes it difficult 
to distinguish between them. We also discuss the case in which the halo density profile is 
fixed independently of the gravity model for comparison between our results and previous 
work. 
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1. Introduction 

Motivated from the discovery of the cosmic accelerated expansion of the univers 
it is becoming important to test the nature of gravity on the cosmological scales. 
This is because the nature of gravity might be deeply rooted to the cosmic acceler- 
ated expansion (e.g., Refs. 02 HI [5]). In general, it is very challenging to construct a 
modified gravity (MG) theory that is consistent with observations. In comparison 
with the varieties of viable dark energy models to explain the accelerated expan- 
sion of the universe, many modified gravity models arc ruled out by cosmological 
observations (e.g., Refs. ES E| El©. 
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For example, one of the mo st p opular modified gravity models is the Dvali- 
Gavadadze-Porrati (DGP) modepl^l. The DGP model is described in the context of 
the brancworld scenario, which consists of a 3+1-dimcnsional branc embedded in a 
5-dimcnsional bulk. This model has an interesting phenomenology and yields two 
branches of the Fricdmann equation, which are a self-accelerating branch (sDGP) 
and a normal branch (nDGP). In the sDGP model, the expansion of the universe 
self- accelerates at late times without a cosmological constant nor spatial curvature, 
while the nDGP model needs to add a stress-energy component with negative pres- 
sure on the brane to be consistent with cosmologic al observa tions. Unfortunately, 
the sDGP model inevitably give rise to a ghost mod e^^ ^^j Moreover, the sDGP 
model is disfavored by the cosmological observation ^ 3 | 14 |15J 

Inspired by the decoupling limit of the DGP model, the galilcon gravity theory 
has been studied as a possible alternative to large distance modification of gravity 
(e.g., Rcfs.[II2[I3[l8l[lIl[Iil[2^ This theory introduces 

a scalar field with the kinetic term with the "wrong" sign and the self- interaction 
term (d<fi) 2 \Dcj), which is invariant under the Galilean shift symmetry d^4> — > d^cfr+b^ 
in the Minkowski space-time, which keeps equation of motion at the second order 
differential equation. Although the Lagrangian no longer satisfies the Galilean shift 
symmetry in a curved spacetime, the field equation of the galileon field remains a 
second order differential equation, and the galilcon field admits the self-accelerating 
solution in a FRW universe without a ghost instability. 

The other notable feature of these modified gravity models is the Vainshtein 
mechanism, which allows the modification of gravity to recover general relativity 
around a high density regional. Thanks to the Vainshtein mechanism, these models 
evade the solar system constraint. It is worth examining whether or not the effect of 
modification of gravity on the nonlinear scales and the structure of the halo of galaxy 
or galaxy cluster might be useful as a clue to test these modified gravity theory. 
Recently, several researchers have investigated such a possibilit y of testing th ese 
modified gravity models on the scales of galaxy or galaxy clusteiElE2l33l34l35] j n 
the present paper, we focus our investigation on the structure of a halo in the sDGP 
model and the galileon model by constructing numerical solutions that correspond 
to the singular isothermal sphere (SIS) and the King model in the limit of Newtonian 
gravity. 

The structure of this paper is as follows. In Sec. [2] we first give a brief review 
of the sDGP model and the galileon model. Then, we give basic formulas for our 
halo models in the sDGP model and the galileon model. In Sec. [3l our numerical 
results of the halo models are demonstrated. In Sec. EJ we investigate the case with 
a different approach for the halo modeling, in which the halo density profile is 
fixed independently of the gravity model for comparison between our results and 
previous work^J Sec. [5] is devoted to summary and conclusions. Throughout this 
paper, we use units in which the speed of light equals unity, and we follow the metric 
convention (—,+,+,+). We use the reduce Planck mass Mpi, which is defined by 
Mpi = 1 /\Z8irG with Newton's gravitational constant G. We adopt the Hubble 
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constant Hq = 100 ft,km/s/Mpc with h = 0.7 and the matter density parameter at 
present Oq = 0.28. 

2. Basic Formulas 

2.1. Modified gravity models 

The Dvali-Gabadadze-Porrati model consists of a spatiall y th ree-dimensional brane 
in a 4+1 dimensional (5D) Minkowski bulk. The action iJ^ 

S = ^ J d 5 x^g~ 5 R 5 + J d A xV=g-(^R + An) , (1) 

where M 5 (-Mpi) is the Planck mass in the 5 (4) dimensional spacctime, g§ (g) and 
i?5 (R) are the determinant and the Ricci scalar of the 5 (4) dimensional metric, 
respectively, and C m stands for the matter Lagrangian. The two Planck masses M 5 
and Mpi can be related via length scale, the crossover scale r c = Mpj/2Mf , which 
must be fine-tuned to be the present-day horizon scales in order to modify gravity 
only at late times. In this model, the Friedmann equation has two branches. One 
corresponds to the self-accelerating solution called sDGP, which contains a ghost- 
like instability. Although the other solution called nDGP does not suffer from a ghost 
instability, the cosmological constant is needed to drive the cosmic acceleration. 
In the present paper, we consider the self-accelerating branch of the DGP model. 
During matter do mination and beyond, the modified Friedmann equation in the 
sDGP model isPl 



H(a) 1-Q Q (1-Op) 2 

-w~^^ + y^ + — i — ' (2) 

where a is the scale factor, and the matter density parameter is related to the 
crossover scale by r c = 1/(1 — Slo)-ffo- Note that the modification of the law of 
gravity in the nDGP model is qualitatively similar to that in the galileon model. 
Therefore, our results in the galileon model can be applied to the nDGP model. 

On the other hand, the galileon model is characterized by a scalar field with 
the self-interaction whose Lagr angian is invariant under the Galilean shift symme- 
try in the Minkowski spacetim c!^^ ! We consider the galileon model in a curved 
spacetime, which is minimally coupled to gravity, with the actiorP^El] 



S = / d^x^g 



M- 



R + K(X) - G{X)U(j> + £„ 



(3) 



where R is the Ricci scalar, X = -g» v V v <j>/2, D<j) = g^V ' v <j>, C m is the 
matter Lagrangian, and K(X) and G(X) is an arbitrary function of X. For sim- 
plicity, we consider the following functions, K (X) = —X and G(X) = (r^/Mp^X, 
where r c is the model parameter. This model admits a late-time de-Sitter attractor 
in a flat FRW universe. The solution along this attractor can remarkably simplify 
the modified Friedmann equation, and the background evolution during the mat- 
ter dominated era can be regarded as the Einstein de-Sitter universe. Besides, this 
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model is not plagued by ghost instability in contrast to the sDGP model. In the 
present paper, we consider the attractor solution. Then, the modified Friedmann 
equation on the attractor can be written as 



(4) 



The parameter r c is related with the cosmological parameters through the relation, 
r c = 1/(54(1 -fio)) 1 /^- 1 - 

Extended models of the galileon model have been propose d I I In Ref. the 
model with G(X) = Mpi(r^/Mp l ) n X n is considered, and it is demonstrated that 
this model reduces to the cosmological constant model for n equal to infinity. Obser- 
vational constraints on this model is also investigated, focusing on the constraints 
from the la supernovae (SN) observations and the cosmic microwave background 
(CMB) distance observation, as well as from the large scale structure of the lu- 
minous red galaxies in the sloan digital sky survey data release 7. The model with 
n = 1 is disfavored by the constraints from the SN and CMB observations, however, 
the model with higher n can be consistent with the observations. In the present pa- 
per, we consider the halo models of the galileon model with n = 1, however, the 
result is almost same as those with higher n at a quantitative level. 



2.2. Perturbation equations in the static approximation 

In this subsection, we summarize basic perturbation equations for gravity and the 
brane bending mode in the sDGP model and the galileon field on small scales, 
assuming the spherically symmetric and static system. For convenience, we choose 
the Newtonian gauge, which is given by 

ds 2 = -(1 + 2*)rft 2 + a 2 (t){l + 2<$>)6 ij dx i dx i , (5) 

where a(t) is the scale factor, but we set a(t) = 1 in the expressions below, for 
simplicity. Following Refs.[l0]and|36]for the DGP model, and Refs.[21]and[29|for the 
galileon model, within the subhorizon scales with the quasi-static approximation, 
we have the following perturbed equations, 

A$ = -AnGp + £Atp, (6) 

$ + * = -aup, (7) 

and 

+ A 2 (vj,«^« - (A^) 2 ) = -4ttGCp, (8) 

where <p(x) denotes the brane bending mode and the perturbation of the galileon 
field defined by <f>(t,x) = <f>(t)(l + <fi(x)), the Laplace operator A represents the 
differentiation with respect to physical space-coordinates, p is the matter density. 
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Here, a, £, £, A 2 and j3 are determined by the background evolution as follows: 

for the sDGP modePES and 

n c a nn 12 m /- Gx4> 2 , 2 G x <t> 



(9) 



f3 = -l + 2G x (4> + 2H<j))-4:TrGG 2 x <j) 4 , (10) 

for the galileon model, where we defined Gx = dG(X) / dX . For the galileon model, 
it is useful to rewrite the combinations and A 2 £ in terms of the matter density 
parameter O m = p m (a)/3Mp l H 2 (a), 

(1 ^m)(2 ^m) -.2/- ( ^ 



which is derived along the attractor solution (see e.g., Ref. [29]) . 

In the spherically symmetric case, Eqs. ©, (JT)) and ([SJ reduce to 

GM(r) d v 

(a + £,) — , (12) 



<ir r 2 tir ' 



<^__r_ 1 . A 8GA 2 CM(r) 
dr ~ 4A 2 I 



(13) 



where we define the enclosed mass M(r) = 4tt J Q r dr'r' 2 p{r'). The influence of </3 IS 
determined by the second term in the square root in Eq. (fT3"]) . which is characterized 
by the so-called Vainshtcin radius, which we defined by 

r v = [8GA 2 CM(rv)] 1/3 . (14) 

On scales smaller than the Vainshtcin radius r « ri/, Eq. (fT3"|) gives dtp/dr = 
— i/GCM(r)/2A 2 r <C GM(r)/r 2 . Therefore, the law of gravity in the Vainshtcin 
limit reduces to Newtonian gravity, 

^ GM(r) 
dr r 2 

On scales larger than the Vainshtcin radius r » ry, Eq. (|13[) gives dip/dr = 
—GC,M(r)/r 2 . Thus, the modification of gravity at large distance becomes 

^ G cff M(r) 
dr r 2 

where G e s = G(l + + is the effective gravitational constant in the linearized 
limit. Geff/G is always less than unity in the sDGP model, while G e ff/G is always 
larger than unity in the galileon model (also in the nDGP model). 
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2.3. Solution of Boltzmann equation 

Next, let us consider the equation that the matter component obeys. Since the 
distribution of the matter is determined through the gravitational potential ty, the 
matter follows the usual collisionless Boltzmann equation. The static solution of 
the Boltzmann equation can be obtained by the Jean's theorem^!. We neglect the 
effect of time-dependent quantities in the modified Poisson equation ([6]) and the field 
equation ©, which might leads to the time-evolution of the effective gravitational 
constant. However, these time-dependent quantities evolve on cosmological time 
scales. Therefore, we can assume this effect can be neglected when considering the 
static solution of a halo. 

In the present paper, we consider the singular isothermal sphere (SIS) model 
and the King model (see e.g., Ref. 137]) . The distribution function is given by 

/w-p^i**"". < 17 > 

for the SIS model, and 

: (is) 

£ < 0, 

for the King model, respectively, where p\ is a constant, £ = — \I> — s 2 /2, s is 
velocity, and a is velocity dispersion. By integrating the distribution function over 
all velocities, we have the formula for the density, 

p = Pl e-*/°\ (19) 

for the SIS model, and 



P = Pi 



U 2 I V 7TC7 2 V 3<T 2 



(20) 



for the King model, respectively, where erf(x) is the error function defined by 
etf(x) = (2/^)J*e- t2 dt. 

The basic equations for a halo are (|T2"j) , (|T3"|) with (fT9"|) for the SIS model, but with 
(f20|) for the King model. These equations are numerically solved in the next section. 
In the case of Newtonian gravity, the solutions arc well knowrJ33. Especially, the 
solution of the SIS model is written in the simple analytic form, p(r) = a 2 /2TrGr 2 . 



3. Numerical results 

3.1. Singular isothermal sphere (SIS) model 

Now, we are ready to solve Eqs. (fT2|) . (fT3|) . and (fT9|) to find the density profile of the 
SIS halo in the sDGP model and the galileon model. Figure Q] shows the density of 
the SIS halo divided by that of Newtonian gravity pgr (r) , as a function of the radius 
r, where PGn( r ) — <J 2 /2nGr 2 . The left panel is the sDGP model, while the right 
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Fig. 1. The left panel is the ratio of the density of the SIS halo in the sDGP model to that of 
Newtonian gravity, as a function of radius r. Here, we adopted the velocity dispersion a = lOkm/s 
(solid curve), <r = lOOkm/s (dashed curve), and cr = lOOOkm/s (dot-dashed curve), respectively. 
The right panel is the same as the left panel, but in the galileon model. 



panel is the galileon model. Each curve corresponds to the velocity dispersion a = 
lOkm/s (solid curve), cr = lOOkm/s (dashed curve), and a = lOOOkm/s (dot-dashed 
curve), respectively. Here, we used the present values of the Hubble parameter and 
the matter density parameter in evaluating (3, A 2 , and £. 

The density ratio, /Cmg/pgRj is almost unity at small radii r in both the models. 
This means that the central region is the same as that of Newtonian gravity. This is 
due to the Vainshtein mechanism. The difference appears as the radius r increases. 
We note that the density ratio finally approaches constant values as r becomes very 
large. Thus, there are two asymptotic regions, i.e., the Vainshtein limit r <C ry and 
the linearized limit r ry. 

In the Vainshtein limit r <C ry, from Eqs. (fT2|) . (|13j) and (jT9j) , we obtain the 
approximate density profile, 




— 7=^r- 

Thus, the effect of modification of gravity is represented by the factor (a + 
OC/V A 2 Cc- On the other hand, in the linearized limit r ^> ry, Eqs. (fTB| and 
(|19[) give the approximate solution 

G 

Pmg(t > ry) ~ —p GR (r). (22) 

Thus, the density is enhanced (decreased) in the sDGP (galileon) model at large 
radii r, reflecting the behavior of G/G e ff = 1/(1 + (a + £)0- 

Let us estimate the Vainshtein radius ry , which can be easily done by substitut- 
ing the relation M = 2a 2 ry /G, extrapolated from Newtonian gravity, into Eq. (|14|) . 
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Then, we have the approximate formula of the Vainshtein radius, 

ry ^ 4y/^a. (23) 

This gives a good approximate expression. However, the exact solution of Eq. (|14[) 
gives a slightly different value from Eq. (|23j) . due to the underestimation (overes- 
timation) of the mass of the halo in the sDGP (galileon) model at r ~ ry. Our 
numerical analysis of Eq. (fT4|) shows that the following formula works, 

r v ~ 2.5 ( f ) ft _1 Mpc, (24) 



lOOkm/s 



for the sDGP model, and 



^^ 4 Hl00kn^J^ lMPC ' ^ 

for the galileon model, respectively. 

We next consider the circular speed of a test particle in a circular orbit at radius 
r, which is given by equating the gravitational force with the centrifugal force 

v 2 (r) = r— . (26) 
dr 

Figure [2] shows the ratio of circular speed in the modified gravity model to that 
in Newtonian gravity, ^mg/^gR' as a function 01 radius r(/i _1 Mpc). Note that the 
circular speed in the SIS halo model satisfies v GR = 2er 2 in Newtonian gravity. The 
left (right) panel of Fig. [5] shows the sDGP (galileon) model. Each curve corresponds 
to the various velocity dispersion a = lOkm/s (solid curve), a = lOOkm/s (dashed 
curve), and a = lOOOkm/s (dot-dashed curve), respectively. 

Inside the Vainshtein radius r < ry, we have v^g ~ v GI{ due to the Vainshtein 
mechanism. However, vj^G becomes to deviate from ^q R as r becomes larger. The 
oscillatory feature appears at r ~ ry. We suppose that the nonlinear structure of 
basic equations plays an important role in the oscillatory behavior around r ~ ry . 
Finally, v^ G approaches Uq R again at large radii r in both models. This behavior 
can be understood by considering the asymptotic solution of the density profile. 
Namely, in the limit r 3> ry, using Eq. (|16|) . we obtain 

4(r»r^ G ^ MG( ^ 4(r), (27) 

with the enclosed mass A/mgM = 47r Jq dr' V 2 pMc( r ') • Since the density of a halo 
in this limit is given by Eq. (|22|) . the effective gravitational constant is canceled 
out in Eq. (|27|) and the circular speed of a test particle in these modified gravity 
models is equal to that of Newtonian gravity. The characteristic behavior at r ~ ry 
appears well outside galaxies or galaxy clusters. Therefore, this effect of modification 
of gravity will be hard to be detected. 
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Fig. 2. The left panel is the ratio of the circular speed of a test particle for the SIS model in the 
sDGP model to that in Newtonian gravity, as a function of radius r(Zi -1 Mpc). Here, we adopted 
the velocity dispersion a = lOkm/s (solid curve), <r = lOOkm/s (dashed curve), and a = lOOOkm/s 
(dot-dashed curve), respectively. The right panel is the same as the left panel, but in the galileon 
model. 




Fig. 3. The density profile of the King model as a function of radius r(h~ Mpc). The density 
is normalized by the critical density at present. The thick curves represent the density profile 
with the central potential \to = —12a 2 , while the thin curves represent the density profile with 
>I/0 = — 3cr 2 . The dashed, solid, and dot-dashed curve correspond to the sDGP model, Newtonian 
gravity, and the galileon model, respectively. Here, the velocity dispersion is fixed as a = lOOkm/s. 



3.2. King model 

The density profile of the King model can be obtained by solving Eqs. (fT2|) . (fT3j) 
and (|20|) . In Fig.[3j we show the density profile of the King model normalized by the 
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critical density at present, as a function of the radius, where the curves correspond 
to the sDGP model (dashed curve), Newtonian gravity (solid curve), and galileon 
model (dot-dashed curve), respectively. We adopted the values of the potential at 
the center, = — 12cr 2 (thick curves) and = — 3<r 2 (thin curves), respectively. 
Here, we set the velocity dispersion a = lOOkm/s. In general, the halo density in the 
King model remains constant at small radii r, and becomes smaller as r increases. 

In the King model, the density of a halo has a cut-off radius r t , called the tidal 
radius, which is defined by ty(rt) = 0. From Eq. (|20j) . we may set p = for r > rt- 
Similar to the case of the SIS halo model, the density is enhanced (decreased) in 
the sDGP (galileon) model in the outer region of halo in comparison with that of 
Newtonian gravity, which explains the behavior of the density near the tidal radius 
in Fig. |3l depending on the models of gravity. 

We define the total mass of a halo by M to t = M(rt)- In the left panel of Fig. HI 
we plot the tidal radius r t as a function of central potential ^>q/<t 2 for the sDGP 
model (dashed curve), Newtonian gravity (solid curve), and the galileon model 
(dot-dashed curve), respectively. The right panel of Fig. |4] plots the total mass 
as a function of the central potential, I& 2 , where each curve is the sDGP model 
(dashed curve), Newtonian gravity (solid curve), and the galileon model (dot-dashed 
curve), respectively. The total mass of the sDGP (galileon) model is larger (smaller) 
than that of Newtonian gravity. The difference of the total mass between the gravity 
models becomes larger as the absolute value of the central potential is larger. 

In the left panel of Fig. 21 we also plot the Vainshtein radius ry. The thin 
dashed curve and the thin dot-dashed curve represent the Vainshtein radius for 
the sDGP model and the galileon model, respectively. In the case r t < ry, the 
density profile at any radius is effectively described by Newtonian gravity, since all 
radii enclosing matter are inside the Vainshtein radius. In the limit of small central 
potential l^ol/V 2 , the tidal radius becomes very small, r t <C ry. In this limit, 
therefore, the difference of the tidal radius between the gravity models becomes 
small. On the other hand, in the limit of large central potential l^ol/c 2 ; the tidal 
radius becomes large, r t ^ ry. In the case r t > ry, the density profile in the outer 
region reflects the modification of gravity. Then, the difference of the tidal radius 
between the gravity models remains in the limit of | Vl/o |/c 2 ^> 1 or r t ry. Let 
us now consider the question whether the gravity model can be distinguishable or 
not. Even for the case r t > ry, the difference in the density profile appears only 
in the outer region of a halo. Thus, the difference between those modified gravity 
theories and general relativity is small, and it will be difficult to be distinguished 
observationally. 

4. Discussion 

In Ref. EU the author investigated the halos in the DGP modified gravity model 
with N-body simulations as well as analytic modeling of a halo similar to our in- 
vestigation. Here, let us discuss the difference between the results in Ref. [31] and 
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Fig. 4. The left panel plots the tidal radius (thick curves) and the Vainshtein radius (thin 
curves) of the King model as a function of 'J'o /<x 2 . Each thick curve is the tidal radius in the sDGP 
model (thick curve), Newtonian gravity (solid curve), and galileon model (dot-dashed curve), 
respectively. The thin dashed (dot-dashed) curve is the Vainshtein radius in the sDGP (galileon) 
model, respectively. The right panel plots the total mass Mtot as a function of ^o/a 2 , for the 
sDGP (dashed curve), Newtonian gravity (solid curve), and galileon model (dot-dashed curve), 
respectively. 



ours. The author in Rcf. (3D has found that the square of the circular speed of a test 
particle is enhanced or decreased in outer region of a halo in proportion to G e g. 
On the contrary, our result in the SIS model is different, which the circular speed 
in the outer region finally approaches the value in the inner region of a halo. The 
difference comes from the treatment of the dark matter component. In our analysis, 
we adopted the explicit formula for the distribution function for the dark matter 
component, which is evidently the solution of the Boltzmann equation. Thus, in our 
approach the same distribution function is assumed for the different gravity model, 
which results in the different density profiles, depending on the gravity models. 

On the other hand, the author of Ref . I3T1 fixed the density profile of a halo inde- 
pendently of the gravity model, which has been suggested from N-body simulations. 
Following this method, the circular speed in a modified gravity model approaches 
that in Newtonian gravity multiplied by the factor G fr/G in the limit of large radii. 
Fig. [5] shows the ratio of the circular speed of a test particle in the sDGP model 
and the galileon model to that in Newtonian gravity with fixing the density profile 
so as to be the SIS model, Pmg i r ) = er 2 /27rGr 2 . In Fig. El we can see the similar 
feature of the circular speed in both the modified gravity models as that expected 
from Rcf. El 

Following Ref. EU we next consider an Navarro- Frenk-White^l (NFW) halo with 
mass M200, defined as the mass contained with a radius -R200 so that the average 
density within i?2oo i s 2 00p, where p is the background density, and we used the 
concentration relation^, c = 9 x (M/[3.2 x 10 12 M Q //i])~ - 13 . In this case, the virial 
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radius is written i?2oo — 1-0 x (il^oo/lO^MQ) 1 / 3 ft. _1 Mpc, and the concentration 
is c ~ 6.0 x (A/ 20 o/10 14 M Q )- 013 (see, e.g., Ref. OS] for details). Fig. [5] shows the 
ratio of the circular speed in the modified gravity models to that in Newtonian 
gravity. Similar feature as that in Fig. [5] can be seen for the NFW halo. In Fig. [5J 
we choose the mass M 200 = 1O 13 M (solid curve), 1O 14 M (dashed curve), 1O 15 M 
(dot-dashed curve), which correspond to clusters. The effect of the modified gravity 
is significant in the outer region of a halo, but it is suppressed in the inner region 
because of the Vainshtein mechanism. For the NFW halo, the Vainshtein radius is 
estimated as 

^"(W*- im *- « 

for the sDGP model, and 

for the galileon model, respectively, from Eq. (|T?1) . The Vainshtein radius is large, 
but the Vainshtein mechanism does not completely hide the effect of the modifica- 
tion of gravity in the cluster. One can read that the circular speed is enhanced or 
decreased in comparison with Newtonian halo, e.g., at 10% level at the radius of 
a few /i _1 Mpc for A/ 2 oo = 10 14 M Q . If a good tracer of circular speed is available, 
there might be a possible chance to test the modified gravity effect in a halo, as 
discussed in Ref. |3U 

5. Summary and conclusions 

In this paper, we have investigated the structure of dark matter halos in the sDGP 
model and the galileon model by using the static and spherically symmetric solutions 
of the collisionlcss Boltzmann equation, which reduce to the SIS model and the King 
model in the limit of Newtonian gravity. We have obtained the solution of a halo in 
these modified gravity theories in a numerical manner. The density of the halo is 
the same as that in Newtonian gravity due to the Vainshtein mechanism well inside 
the Vainshtein radius r <C ry. We have found that the density in the sDGP (the 
galileon) model becomes larger (smaller) in the outer region of a halo, in comparison 
with that in Newtonian gravity, which comes from the suppression (enhancement) 
of the effective gravitational constant in the sDGP (galileon) model, respectively. In 
the SIS model, the density of the halo at large radii r ^ ry approaches the value of 
Newtonian gravity multiplied by a factor G /G c ff- The circular speed is also modified 
around r ^ ry. However, the Vainshtein radius is ry ~ 1 ~ lOMpc for a typical 
galaxy and ry ~ 10 ~ lOOMpc for a typical galaxy cluster. Therefore, it will be 
difficult to distinguish between these modified gravity theories by a measurement 
of a halo because the Vainshtein radius is large enough. 

However, the above results rely on the assumption that the distribution function 
of the dark matter follows the simple explicit formula. For a comparison with the 
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Fig. 5. The ratio of the circular speed of a test particle in the galilcon model (thick curves) 
and the sDGP model (thin curves) to that in Newtonian gravity for the SIS model, as a function 
of radius r(h~ Mpc), when the density profile of a halo is fixed so as to be the SIS model, 
Pmg i r ) = o- 2 /2irGr 2 . Here, we adopted the velocity dispersion a = lOkm/s (solid curve), o = 
lOOkm/s (dashed curve), and a = lOOOkm/s (dot-dashed curve), respectively. For large r, this 
ratio approaches G c ff/G = 1 + £C — 1-94 for the galileon model, and G a g/G = 1 + 1/(3/3) ~ 0.715 
for the sDGP model, respectively. 




Fig. 6. The same as Fig. [5] but for the NFW halo. Here, we adopted the virial mass M200 = 
10 13 M Q (solid curve), 10 14 Mq (dashed curve), and 10 15 Mq (dot-dashed curve), respectively. The 
asymptotic value of G e ff /G at large r is same as that in Fig. [5] 

previous approach in Ref. 131 1 we investigated the effect of the modification of gravity 
with fixing the halo density profile independently of the gravity model. In this case, 
the circular speed is enhanced or decreased in the outer region of a halo depending 
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on the gravity model. In the case adopting the NFW profile as a model of a cluster 
halo, the Vainshtein radius is large. However, the Vainshtein mechanism does not 
completely hide the effect of the modification of gravity in the cluster. This might 
provide a possible chance that precise measurement of halo could be a probe of the 
modified gravity, as discussed in Ref. [3T1 
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